Long-range correlation of thermal radiation 
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A general theory is presented for the spatial correlations in the intensity of the radiation emitted 
by a random medium in thermal equilibrium. We find that a non-zero correlation persists over 
distances large compared to the transverse coherence length of the thermal radiation. This long- 
range correlation vanishes in the limit of an ideal black body. We analyze two types of systems (a 
disordered waveguide and an optical cavity with chaotic scattering) where it should be observable. 
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The Hanbury Brown-Twiss effect is the existence of 
spatial correlations in the intensity of thermal radiation 
by a distant source. It was originally proposed as an 
intensity-interferometric method to measure the angular 
opening of a star |jj , far less susceptible to atmospheric 
distortion than amplitude- interferometric methods |Q. 
Two photodetectors at equal distance r from a source (di- 
ameter a) will measure a correlated current if their sep- 
aration d is smaller than the transverse coherence length 
d c ~ Xr/a of the radiation from the source at wavelength 
A. The correlation function decays with increasing d in 
an oscillatory way, with amplitude oc (d c /d) 3 

The textbook results assume that the source of the 
thermal radiation is a black body, meaning that at each 
frequency any incident radiation is either fully absorbed 
or fully reflected. In a realistic system there will be a 
frequency range where only partial absorption occurs. 
The purpose of this paper is to show that in general for 
thermal radiation the correlation function does not de- 
cay completely to zero, but to a non-zero d-independcnt 
background value. This long-range correlation is smaller 
than the short-range correlation by a factor (A/a) 2 , and 
becomes dominant for d > r{\/ a) 1 ^ . It contains infor- 
mation on deviations of the thermal radiation from the 
black-body limit. 

The new information contained in the long-range cor- 
relation is most easily described when the source is em- 





FIG. 1. Schematic diagram of a source (length L, diameter 
a) radiating into an iV-mode waveguide that is open at both 
ends. The radiation leaving the waveguide at one end is de- 
tected by two photodetectors at a distance r from the source 
and separated by a distance d. The photocathodes have an 
area below the transverse coherence area d 2 c ~ r 2 /N. We find 
that the photocurrents are correlated even if the two detectors 
are separated by more than do- 



bedded in a waveguide (see Fig. [I]). The waveguide has 
length L, cross-sectional area A ~ a 2 , and supports 
TV = 2ttA/X 2 propagating modes at frequency ui, count- 
ing both polarizations. In the far-field, and near nor- 
mal incidence, each mode corresponds to a transverse 
coherence area (r\) 2 /A = d 2 . The source is in ther- 
mal equilibrium at temperature T . The radiation emit- 
ted through the left end of the waveguide is incident on 
a pair of photodetectors, one detecting the photocur- 
rent Ik{t) in mode k, the other detecting J/(i). Each 
photocathode has an area equal to the coherence area 
or smaller. The photocount rik = rik + Srik (number 
of photons counted in a time t) and the photocurrent 
Ik = dnk/dt = Ik + 5 Ik fluctuate around their time- 
averaged values rik and Ik = u k jt. We seek the correla- 
tion function 
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C k i = / 5I k {t + T)5h{t) dr = lim -5n k (t)6m(t) . 

(1) 

The overline indicates an average over many measure- 
ments on the same sample. 

The advantage of embedding the source in a waveguide 
is that we can characterize it by a finite-dimensional scat- 
tering matrix S(u>), consisting of four blocks of dimension 
N x N, 



S 



r t 
t' r' 



(2) 
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A mode I incident from the left is reflected into mode k 
with amplitude and transmitted with amplitude t' kl . 
Similarly, r' kl and tki are the reflection and transmission 
amplitudes for a mode I incident from the right. Reci- 
procity relates these amplitudes by rjy = rik, r 'ki 
and tki = t' lk . 

It has been shown recently by one of the authors [Q, 
using the method of "input-output relations" how 
the photocount distribution can be expressed in terms of 
the scattering matrix. The expressions in Ref. B are for 
a single multi-mode photodetector. The corresponding 
formulas for two single-mode photodetectors are 



C kl =ockOLi\ \(QQ r )ki{u)\ z [f{io,T)] 
Jo 
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h=oc k \ (QQ^) kk (ui)f(ui,T)^- 

Z7T 



(3) 



where a k is the detector efficiency (the fraction of the 
photocurrent in mode k that is detected) and / is the 
Bose-Einstein function 



/(w,T) = [exp(hu/k B T) - If 1 



(4) 



The N x N matrix Q is related to the reflection and 
transmission matrices by 



QQ f = 1 - rr f - tt 1 



(5) 



The integral over u> extends over a range Q c set by the 
absorption line width, centered at ujq. Typically, Q c <C 
lu , so we can neglect the frequency dependence of N and 
/. The matrix Q(u) for a random medium fluctuates on 
a scale uj c much smaller than Q c . The integration over 
u) then averages out the fluctuations, so that we may 
replace the integrand by its ensemble average, indicated 
by (...), 



Cu = ctkaif 



{\{QCp)u{u)\ 2 )—+8 kl Ik- (6) 
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We evaluate the ensemble average using results from 
random-matrix theory Q . For a medium with randomly 
placed scatterers, the "equivalent channel approxima- 
tion" H has proven to be reliable. According to this 
approximation, all N modes are statistically equivalent. 
As a consequence, for any k ^ I one has 



/V 



(tr(QQ 



t\2\ 



= Nj2((QQt) kj (QQl) jk ) 

= N(N - 1)(\(QQ%\ 2 ) + N((QQi)t k ) . (7) 
The average of (QQ^) 2 k factorizes in the large- N limit || , 

((QQ^)l k ) = {(QQ% k ) 2 [l + 0(N- 1 )] 

= 7V- 2 (trQQt) 2 . (8) 

Combination of Eqs. (^) and (^) gives us 

(\(QQ%\ 2 ) = N- 2 (tv(QQ^) 2 ) 

- JV _3 (tr QQ + ) a + 0(N- 2 ) . (9) 

The eigenvalues oi, 02, . . . , <tn of the matrix rr^ + tt^ 
are the "scattering strengths" of the random medium. 
We denote by ~&p = N^ 1 ^ n a? the p-th. spectral moment 
of the scattering strengths. According to Eqs. (||), (||) 
and (^|), the cross-correlator C k i (k ^ I) then takes the 
form of a variance, 



Cu — 



N 



did 

2^ 



This is our basic result for the long-range correlation an- 
nounced in the introduction. The new information con- 
tained in the cross-correlator is the variance of the scat- 
tering strengths. The auto-correlator, in contrast, de- 
pends entirely on the first spectral moment, 



C kk = a 2 J 2 ^(l-wf^+h 



(l-a 



2^ 
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where we have used Eq. (g). 

The long-range correlation C k i of two photodetectors 
separated by more than a coherence length is an order N 
smaller than the short-range correlation Ckk — Ik of two 
photodetectors separated by less than a coherence length. 
(The full value Ckk is measured in a single-detector ex- 
periment.) The long-range correlation vanishes if all N 
scattering strengths are the same, as they would be for an 
idealized "step-function model" of a black body (a n = 
for \u> — luq\ < f2 c and a n = 1 otherwise). A random, 
partially absorbing medium, in contrast, has a broad dis- 
tribution of scattering strengths S, hence a substantial 
long-range correlation of the photocurrent. 

As first example, we compute the correlation for a 
weakly absorbing, strongly disordered medium. The mo- 
ments of rr^ and ttf appearing in Eqs. ( |Io|) and ( |TT| ) have 
been calculated by Brouwer [jlO| as a function of the num- 
ber of modes N, the sample length L, the mean free path 
I, and the absorption length £ = \/Dr a (r a is the absorp- 
tion time and D — cl/3 the diffusion constant). It is 
assumed that 1/N <C 2/£ <C 1, but the ratio L/£ = s is 
arbitrary. The result is 



n- ) - (a) 2 = — coth 3 s - — - + — — 
3£ V sinhs smh s 

s coth s — 1 s 



sinh 

/, 4/ , S 

(1 — a) = — tanh - 
V 1 3£ 2 



sinh 
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To compute the correlators (|10|) and flll| ) it remains 
to carry out the integrations over u>. The frequency de- 
pendence is governed by the imaginary part of the di- 
electric function e"(u>), for which take the Lorentzian 
e"(u) = e' '[l + (w - ljo) 2 /^}- 1 . Since r a = l/u> e", the 
corresponding w-dependence of £ and s is £/£o = sq/s — 
[1 + (lj — luq) 2 /SI 2 ] 1 / 2 , with £ and so the values of £ and 
s at u> = ujq. Results are plotted in Fig. ||. In the limit 
L/£o — > of a thin sample, we have 

Cu = (ifWat/NZo) (L/£ ) 3 , (13a) 
Ckk = (lfa k /Zof (L/^) 2 + Ik , (13b) 

Ik = (iM/Co) • (13c) 



2 




FIG. 2. Long-range correlation Cki (solid line), in units of 
Q. c lf 2 akCii/N^o, and short-range correlation Ckk — Ik (dashed 
line), in units of fi c (Z/afc/£o) 2 , of the radiation emitted from 
a disordered waveguide (inset). A Lorentzian frequency de- 
pendance is assumed for the dielectric function, with width 
£l c and absorption length £o a t the center of the absorption 
line. The mean free path I is assumed to be <C £o. The 
short-range correlation saturates in the limit L/£o — > oo, 
while the long-range correlation keeps increasing oc lnL/£o- 



In the opposite limit L/^o ^ oo of a thick sample, the 
cross-correlator Cki and the mean current Ik both diverge 
logarithmically oc ln£/£o- The ratio Cki/ (Ikh) 1 ^ 2 tends 
to (l/2N)fy/ctk(xi in the large- L limit, and the short- 
range correlation C k k — Ik to ^ c (U a k/^o) 2 , remaining 
larger than the long-range correlation because the limit 
TV — > oo has to be taken before L — ► oo. 

Our second example is an optical cavity filled with an 
absorbing random medium (see Fig. ||a, inset). The ra- 
diation leaves the cavity through a waveguide support- 
ing N modes. The general formula (|§|) applies with 
QQ^ = 1 — rr* (since there is no transmission). The 
scattering strengths a\, . . . , ctjv in this case are eigen- 
values of rr^ . Their distribution is known in the large- N 
limit Q as a function of the dimensionless absorption 
rate 7 = 2?r /Nr a Auj, with Auj the spacing of the cavity 
modes near frequency luq. (The quantity 7 is the ratio 
of the mean dwell time in the cavity without absorption 
and the absorption time.) The moments (a) and (a 2 ) 
can then be computed by numerical integration. Results 
are shown in Fig. again for a Lorentzian frequency de- 
pendence of e"{uj). Unlike in the first example, we are 
now not restricted to weak absorption but can let the 
absorption rate 70 at the central frequency ujq become 
arbitrarily large. For weak absorption, 70 -C 1, we have 



Cki = \n c (fa k ai/N)y$ , 
C'kk = ^c(/a fc 7o) 2 +Ifc, J k 



iri c /cv fc 7o 



(14a) 
(14b) 



Cki/ilkhYi 2 tends to 0.062f(a k aiY/ 2 /N, and the ra- 
tio (C'kk - h) I h to \ja k (see Fig. |b). The long-range 
correlation does not vanish as 70 — > 00, because there re- 
mains a tail of frequencies with moderate absorption and 
thus a wide distribution of scattering strengths, even if 
the system behaves like an ideal black body for frequen- 
cies near ojq. 

In summary, we have shown that the thermal radia- 
tion emitted by random media contains long-range spa- 
tial correlations in the intensity. The long-range corre- 
lation has information on the spectral variation of the 
scattering strengths that is not accessible from the lumi- 
nosity. We have analyzed two types of systems in detail, 
providing specific predictions that we hope will motivate 
an experimental search for the long-range correlation. 

This work was supported by the Dutch Science Foun- 
dation NWO/FOM. 
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For strong absorption, 70 S> 1, all three quantities C k i, 
Ckk and Ik diverge oc ^/to (see Fig. ||a). The ratio 
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FIG. 3. Correlators of the radiation emitted from a disor- 
dered optical cavity (inset) as a function of the absorption 
rate 70 at the center of the absorption line with Lorentzian 
profile. (The absorption rate is normalized to the mean dwell 
time.) (a) Long-range correlation Cm (solid line), in units of 
fl c f 2 akCti/N, and short-range correlation Ckk — Ik (dashed 
line), in units of Q, c f 2 a\. (b) Same correlators, but now nor- 
malized by the mean photocurrent. (The left axis is in units 
of fy/oikCti/N, the right axis in units of /a&.) The long-range 
correlation persists in the limit 70 — > 00 because of partial 
absorption in the tails of the absorption line. 
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